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Abstract 

Electron transport in bilayer graphene is studied by using a first principles analysis and the Monte 
Carlo simulation under conditions relevant to potential applications. While the intrinsic properties 
are found to be much less desirable in bilayer than in monolayer graphene, with significantly 
reduced mobilities and saturation velocities, the calculation also reveals the dominant influence 
of extrinsic factors such as the substrate and impurities. Accordingly, the difference between 
two graphene forms are more muted in realistic settings although the velocity-field characteristics 
remain substantially lower in the bilayer. When bilayer graphene is subject to an interlayer bias, 
the resulting changes in the energy dispersion lead to stronger electron scattering at the bottom of 
the conduction band. The mobility decreases significantly with the size of the generated bandgap, 
whereas the saturation velocity remains largely unaffected. 

PACS numbers: 72.80.Vp,72.10.Di,72.20.Ht,73.50.Dn 
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I. INTRODUCTION 



Graphene has received much attention in the last few years due to its unique properties. 
In addition to the significant interests in fundamental physics, stemming in part from the 
relativistic-like behavior of the massless charge particles around the Dirac cone, this mate- 
rial is very attractive in many applications, particularly in high speed devices.- - - However, 
the gapless electron spectrum of monolayer graphene (MLG) makes it difficult to turn off 
the electrical current due to tunneling. Bilayer graphene (BLG), on the other hand, can 
provide a finite band gap up to hundreds of meV, when the inversion symmetry between top 
and bottom layers is broken by an applied perpendicular electric field." A current on/off 
ratio of about 100 was observed at room temperature,- offering a much needed control for 
nonlinear functionality. Unfortunately experimental studies have also indicated that the 
typical mobility of charge carriers in BLG may be substantially smaller than in MLG.— A 
recent work based on the first principles calculations has suggested that this discrepancy 
may start with the intrinsic transport properties, which results from substantial differences 
in electron-phonon coupling in these two materials.— Additionally, it is widely accepted that 
extrinsic factors such as charged impurities, disorder, and surface polar phonons (SPPs) can 
significantly alter carrier transport in graphene on a substrate, a most commonly used con- 
figuration.— _ — Despite extensive research efforts, a comprehensive understanding of electron 
transport properties in BLG is still a work in progress. 

The purpose of the present paper is to address this issue theoretically by taking advantage 
of a first principles analysis and the Monte Carlo simulation. Specifically, the impacts of 
substrate conditions and perpendicular electric fields are examined. The calculation results 
indicate that graphene in the bilayer form loses much of its advantage over conventional 
semiconductors in the low field transport, particularly when the band structure is modified 
to induce non-zero energy gap. The saturation drift velocity, on the other hand, can re- 
main relatively high. Due to the screening, electrons in BLG appear less susceptible to the 
interactions with remote Coulomb sources, such as SPPs and impurities on the substrates, 
compared to the monolayer counterparts. Below, we begin with a brief overview of the 
models used to estimate the relevant scattering rates. 
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II. RELEVANT SCATTERING MECHANISMS 



The strength of the electron-phonon coupling is estimated from the first principles by 
using the density functional theory.—^ A comparative analysis reveals several qualitative 
differences in the intrinsic scattering of MLG and BLG. For one, MLG has six phonon 
branches with two carbon atoms in a unit cell, whereas these numbers double in BLG. 
Then, BLG may need to consider both the intraband and interband transitions due to 
the presence of a close second conduction band 7r% in addition to the lowest n{ states. In 
addition, the optical phonons in BLG appear to be a relatively weak source of interaction 
unlike in MLG. Ultimately, the intrinsic scattering rate in BLG is dominated by the long 
wave acoustic phonons (intravalley scattering). Figure [T] illustrates this general trend; only 
the dominant branches are shown for clarity of presentation. The nomenclature for the 
phonon modes can be found in Refs. and Il5. 



In the presence of a polar substrate, the graphene electron interaction with SPPs can 
play a significant role as the earlier studies in MLG have illustrated comprehensively.— 1 ^ 1 ^ 
A similar treatment can be extended to BLG. By assuming that the electrons are equally 
distributed in the two layers of BLG, we can derive the corresponding scattering rate as 
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where q = |kj — kj| is the magnitude of the SPP momentum, Ef (Ei) is the final (initial) 
electron energy, d is the distance between the first layer and the substrate (0.4 nm), c 
is the interlayer distance (0.34 nm), us is the SPP energy, n q is the phonon occupation 
number, J 72 = ^j^( K o} +1 — is the square of Frohlich coupling constant, and e(q) is the 
dielectric function. The term | g^ s (q)\ 2 = | [1 + ss' cos Ok cos «k+q + ss' sin at sin ctk+q cos 26] 
comes from the overlap of the electron wave functions of the initial and final states with the 
scattering angle s and s' are the band indices whose product is +1 for the intraband 
(for example, vr^-vr*) and —1 for the interband (7r^-7r2) transitions. 1 - For intrinsic BLG, 
at = vr/2 for an arbitrary k. When an interlayer bias of u is applied, it is modified to satisfy 
tan«k = h 2 k 2 /m*u, where m* is the effective mass of unbiased (or intrinsic) BLG at the K 
or K' point.— For MLG, |g£' s (q)\ 2 = |[1 + ss' cos 8}. The static dielectric function e(q) in 
BLG at room temperature can be estimated by using the random phase approximation. In 



addition, T contains the dependence on the high (low) frequency dielectric constant of the 
substrate The specific values for the relevant substrate parameters can be found 



in Refs. 
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As for the remote impurity scattering, the charged impurities are 
assumed to be located on the surface of the substrate, in which case the scattering rate is 
given by 
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where k = (k° s + l)/2 is the background dielectric constant and rii is the impurity concen- 
tration. The other parameters are the same as defined in Eq. (Op). 

The dielectric function in BLG can be expressed as^^i e(q) = 1 — v q Il(q) with the bare 
Coulomb interaction v q = e 2 /2e q and the electron- hole propagator 

m = 2j2\ 9 s /(q)\ 2 f Jr~_ f ; s , (3) 
S)S ',k k +q k 

where is the electron distribution function in band s and the factor of 2 takes into account 
the spin degeneracy. Figure [2] shows the numerically obtained propagators for MLG and 
BLG with the graphene electron density of n = 5 x 10 11 cm -2 at 300 K. In this calculation 
[i.e., 11(g)], electron occupation in the 7r| states is ignored for its negligible contribution. 
Compared to MLG, the screening in BLG appears to be much stronger due mainly to the 
large density of state at the bottom of the conduction band (ttI).— Additionally, the impact 
of the interlayer potential on the electron screening in BLG can be substantial and is a 
strong function of temperature. 



III. ELECTRON TRANSPORT IN BLG VS. MLG 

A full-band ensemble Monte Carlo calculation is used for evaluating electron transport 
characteristics self-consistently. The model takes into account the complete electron and 
phonon spectra in the first Brillouin zone. Specifically, both the graphene phonon disper- 
sion and its interaction with electrons are obtained from the first principles calculations as 
discussed above,— >^ whereas analytical expressions from the tight binding approximation 
are used for the electronic energy bands in MLG and BLG (with the nearest-neighbor hop- 
ping energy 7 = 3.3 eV and and the interlayer hopping energy 71 = 0.4 eV).— The effect 
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of degeneracy in the electronic system is taken into account by the rejection technique, af- 
ter the final state selection. Electron scattering with the SPPs and remote impurities are 
also considered in the calculation as described above whenever necessary. Throughout the 
calculation, we assume n = 5 x 10 11 cm -2 and T = 300 K. 

Figure [3] shows the electron drift velocity as a function of the electric field in MLG and 
BLG (with no interlayer bias). It provides a comparison of all examined cases: namely, 
intrinsic graphene and graphene on two different substrates, SiC>2 and BN, for which a 
charged impurity density of rii = 5 x 10 11 cm -2 is considered along with the SPPs. As 
illustrated, intrinsic MLG (i.e., no substrate) shows remarkable transport properties,—* 2 ^ 
with the mobility and the saturation velocity as high as 1.0 x 10 6 cm 2 /Vs and 4.3 x 10 7 cm/s. 
An analogous calculation for BLG gives a substantially lower mobility of 1.2 x 10 5 cm 2 /Vs 
and the saturation velocity of 1.8 x 10 7 cm/s. These results for the intrinsic drift velocity 
are consistent with the scattering rates in Fig. [lj they demonstrate how acoustic and optical 
phonons affect the charge carriers in MLG and BLG at various electric fields. 

The lower mobility in BLG can be readily explained by the larger acoustic phonon scat- 
tering rates, as well as the lower electron speed near the bottom of the conduction band 
where all electrons tend to congregate at low electric fields. To understand the smaller 
saturation velocity in BLG, on the other hand, it is convenient to examine the distribution 
function at high electric field plotted in Fig. HI As the electrons gain energy in the applied 
field, the distribution function shifts in the k-space along the direction of the drift. At 
the same time, increased scattering with the long-wave acoustic phonons leads to a further 
broadening of the electron distribution. The result of this interplay between the displace- 
ment and the broadening is the saturation of the drift velocity in intrinsic BLG, where the 
acoustic phonons dominate the scattering. In MLG, however, the velocity curve appears to 
demonstrate another pattern, which points to the presence of a competing intrinsic scat- 
tered Indeed, as we discussed earlier, the inelastic scattering via optical phonons is strong 
unlike in BLG, providing energy relaxation for hot electrons. Consequently, the distribution 
in MLG is prevented from shifting further towards higher energies, which in turn effectively 
curtails the momentum relaxing interactions and results in a higher saturation velocity. 

A similar consideration applies when an additional source of optical phonon scattering 
comes into play. That is, the saturation velocity may be actually enhanced when graphene 
electrons are subject to the SPP scattering, as it provides another path for hot electron 
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energy relaxation. Figure |3] clearly demonstrates this effect. In MLG, the SPP scattering 
can increase the saturation velocity up to 6.5 x 10 7 cm/s if SiC>2 or BN is used for the 
substrate material. This trend of the positive impact of the substrate on the saturation 
velocity was confirmed in recent studies.— In BLG, the saturation velocity can reach 
2.9 x 10 7 cm/s, which is about as 1.5 times large as the intrinsic value. Even though the 
screening in BLG is larger, leading to smaller SPP scattering rates, the enhancement of 
drift velocity is still prominent. This is due partly to the fact that the competing relaxation 
mechanism (i.e., intrinsic optical phonon scattering) is weak in BLG as discussed earlier. 

When the charged impurity scattering is taken in account, the electron drift velocities in 
MLG and BLG are substantially reduced at low electric fields (see Fig. [3]). If BN is used as 
the substrate, the low field mobility is 1.9 x 10 4 cm 2 /Vs in MLG, and 1.2 x 10 4 cm 2 /Vs in 
BLG. The drift velocity at the electric field of 20 kV/cm in MLG decreases to 4.8 x 10 7 cm/s, 
while in BLG it stays nearly the same as the case without impurity scattering, 2.8 x 10 7 cm/s. 
Clearly, the impact of ionized impurity scattering is much less severe in BLG. This is because 
the BLG electrons on average are more separated from the impurities (i.e., the surface of the 
substrate) and experience stronger screening. The calculated result in MLG on Si02 shows 
a good agreement with the available experimental data;— the discrepancy may be attributed 
to the presence of additional sources of interaction such as neutral scatters on the substrate, 
ripples and other defects in graphene crystal lattice.- 

IV. TRANSPORT IN BLG WITH INTERLAYER BIAS 

BLG with an interlayer bias offers the advantage of a tunable bandgap with potential 
applications to transistors, tunable photo-detectors and lasers.->2 The band structure in this 
case changes to, 5 

El = | + £ + nW F tf ± i/| + ^( 7 f + ^), (4) 

where u is the difference between on-site energies in the two layers, vf = ( a/3/2) Gryo/fr ~ 
1 x 10 s cm/s is the Fermi velocity (a = 0.246 nm), and the minus and plus signs correspond 
the 7rJ and conduction bands, respectively. However, this ability to tune the band gap 
comes at the expense of the material's intrinsic transport properties. As the bandgap opens, 
the bottom of the lowest conduction band changes its shape from a hyperbola to a so-called 
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"Mexican hat" and the density of states exhibits a Van Hove singularity as shown in Fig. [51 
The inserted figure shows the band structures for different bandgaps. Consequently, the 
increased density of states leads to stronger quasi-elastic electron interaction with the long- 
wavelength acoustic phonons at low electron energies. Figure M provides the dependence of 
the mobility on the size of the bandgap. When the bandgap is 0.24 eV, the mobility drops 
as low as to 1.2 x 10 4 cm 2 /Vs even without any external scattering mechanisms, which is 
one order of magnitude smaller than in the unbiased case. The reduction in the mobility 
becomes even more pronounced when the ionized impurity scattering is taken into account. 
The Van Hove singularity also enhances the impact of electron coupling with the remote 
impurities. 

On the other hand, it appears that the velocities at high fields are mostly unaffected by 
the gap. Figure [7J shows the drift velocity for u = 0.1 V, 0.2 V, and 0.3 V, which correspond 
to the bandgap of 0.1 eV, 0.18 eV, and 0.24 eV, respectively. As the velocity saturation 
is associated with hot electrons, the changes of the electron spectrum on the bottom of 
the conduction band have little to no effect on the charge carriers located much higher 
in energies. Similarly, the impact of the SPP scattering on the drift velocity, as it is felt 
primarily via the high energy electrons, is not affected by the gap in the electron energy 
spectrum. Figure [7] demonstrates this point as no appreciable difference is observed in the 
SPP contribution for three different choices of u. When the ionized impurity scattering is 
included, the velocity saturation is progressively pushed to a higher field due to the reduction 
in the mobility (i.e., the slope) that was discussed above. 

It should be noted that the mobilities obtained in the presented Monte Carlo simulation 
may have a limited accuracy, for lack of an available adequate model for the electron-electron 
scattering. It might be especially important in case of biased BLG, where the divergence 
of the density of states results in difficulties of reproducing electron thermal equilibrium in 
Monte Carlo simulation. 

V. SUMMARY 

Electron transport properties of BLG are studied under realistic conditions in the presence 
of the SPPs and charged impurities. Overall, BLG has a lower mobility and saturation 
velocity than MLG, due to the stronger acoustic phonon scattering, weaker optical phonon 
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scattering, and nonlinear dispersion at the bottom of the conduction band. It is also shown 
that SPPs can improve the saturation velocity in BLG by effectively dissipating the electron 
energy The impurity scattering has a strong effect of decreasing the drift velocities in both 
MLG and BLG. However, BLG appears more resistant to impurity scattering than MLG, 
due to a stronger screening and larger effective distance between electrons and the impurities. 
In BLG with a interlayer bias, the changes in the band structure result in drastically reduced 
mobilities, especially in the presence of charged impurities. This may limit the chances of 
utilizing BLG in device applications. 
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FIG. 1: (Color online) Intrinsic electron scattering rates in MLG and BLG calculated from the first 
principles for the dominant (a) acoustic and (b) optical phonons . lllls The sudden increases shown 
in (b) are due to the onset of optical phonon emission. They correspond to the phonon frequencies 
at the points of high symmetry in the first Brillouin zone: uj-p ~ 200 meV and ujk ~ 160 meV. 
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FIG. 2: (Color online) Electron-hole propagator Tl(q) normalized to Nq (= 2m* /irh 2 ) in MLG and 
BLG with the graphene electron density of 5 x 10 11 cm -2 . For BLG, the calculation considers 
different interlayer biases with the induced energy gap E g of eV (no bias), 0.10 eV, 0.18 eV, and 
0.24 eV, respectively. 
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FIG. 3: (Color online) Electron drift velocity versus electric field in (a) MLG and (b) BLG, with 
different substrate conditions: intrinsic/no substrate (circle), Si02 (triangle), Si02 with impurities 
(reverse triangle), BN (square), and BN with impurities (diamond). The electron density is 5 x 
10 11 cm -2 at 300 K. The impurities on the surface of the substrate (d=0.4 nm) have the density 
5 x 10 11 cm" 2 . 



13 




k x (10 7 cm" 1 ) 

FIG. 4: (Color online) Cross sectional view (£^=0) of electron distribution functions in (a) MLG 
and (b) BLG at 20 kV/cm, with different substrate conditions: intrinsic/no substrate (short dashed 
line), Si02 (long dashed line), Si02 with impurities (solid line), BN (dashed-dotted line), and BN 
with impurities (dotted line). The conditions are the same as specified in Fig. [3]. 
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FIG. 5: (Color online) Density of states in BLG for different interlayer biases with the induced 
energy gap E g of eV (solid line), 0.10 eV (dashed line), 0.18 eV (dotted line), and 0.24 eV 
(dashed-dotted line), respectively. In the insert, the corresponding band structures are plotted. 
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FIG. 6: (Color online) Electron mobility versus bias-induced bandgap in BLG at 300 K, with 
different substrate conditions. The electron density is 5 x 10 11 cm -2 . The impurities on the surface 
of the substrate (d=0.4 nm) have the density 5 x 10 11 cm~ 2 . 
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FIG. 7: (Color online) Electron drift velocity versus electric field in BLG with a bias-induced gap 
of Eg of (a) 0.1 eV, (b) 0.18 eV, and (c) 0.24 eV, under different substrate conditions: intrinsic/no 
substrate (circle), Si02 (triangle), Si02 with impurities (reverse triangle), BN (square), and BN 
with impurities (diamond). The conditions are the same as specified in Fig. [6J 
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